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Abstract 

We have measured the surface impedance of thick superconductors in the mixed 
state over a broad 2 kHz - 20 MHz frequency range. The depinning cross-over is 
observed; but it is much broader than expected from classical theories of pinning. 
A striking result is the existence of size effects which invalidate the common inter- 
pretation of the low-frequency surface inductance in terms of a single penetration 
depth. Instead, a two-mode description of vortex dynamics, assuming free vortex 
flow in the bulk and surface pinning, accounts quantitatively for the spectrum of 
the complex apparent penetration depth. 

Key words: vortices, superconductor, surface impedance, pinning, depinning 
transition. 



1 Introduction 



The presence of an array of quantized vortices (ip = lines) in a type-II super- 
conductor entails a frictional force in the Euler equation for the supercurrent. 
For an isotropic material it takes the form [1] 

E H - 1 = (3u: x [( J s + curl e) x u] + (3'uj x ( J s + curl e) = w x v L , (1) 

6 Ob 

where V s , J s = —n s eV s and n s are coarse-grained mean values of the superfmid 
velocity, supercurrent, and density of superelectrons; m and e are the electronic 
mass and charge. The electric field E is directly related to the vortex field uj = 
riyip^v and the line velocity v L through the generalized Josephson equation. 
The vector uj incorporates both density riy and orientation v of the vortex 
lines. 
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The macroscopic London equation, [2,3] 
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us = B curlV s = B + curl(/i A 2 J s ), (2) 

e 

shows that vortex lines may deviate from magnetic field lines (u> ^ B) in 
the presence of a supercurrent. In Eq. (2) A is a field-dependent penetration 
depth. 

Here we consider a slab in a perpendicular applied field B . The vortex density 
ny — -Bo/Vo is nearly constant, but v may be quite different from B/B; such 
bending effects are particularly important near the sample surface. 

Eq. (I) is the superconducting analog of the acceleration equation for rotat- 
ing superfluids (Eq.(24) in Ref. [4]) and has been written in its most general 
form such as deduced from conservation laws. The line tension e = ev is the 
thermodynamic conjugate variable of uj. Kinetic coefficients (3 and f} 1 describ- 
ing the longitudinal (resistive) and transverse (reactive) components of the 
friction force are analogs of the mutual friction parameters B and B' in ro- 
tating superfluids [5,6]. Allowing for a normal-current contribution J n = aE 
(a = <7||,<7' = <7j_) in the total current J = J n + J s , (3 is related to the longi- 
tudinal flux-flow conductivity Of = er+1/ '/3u. The Hall-angle is related to f3' 
through ttm6 H = (a' + l/[3'uj)/a i . 

In conventional dirty superconductors dissipation (f3) arises from normal- 
current back-flow and relaxation phenomena at the vortex core, while Hall 
angles are small. Of special interest for future experiments will be the case 
of clean superconductors at low temperatures where Of ~ 1//3uj and Oh are 
expected to be governed by the discreteness of core levels in relation with the 
symmetry of the order parameter [7,8]. 

Vortex dynamics is best probed by measuring the surface impedance Z = 
R — iX in the linear regime, which is conveniently expressed in terms of an 
effective complex penetration depth 

A perfect sample, completely free from defects, should behave like a contin- 
uous resistive medium with Z = (I + i)5f/2, where 5f = (2//ioO"f^) 1 ^ 2 is the 
usual skin-depth. In real samples, this ideal response has been observed in 
conventional superconductors at microwave frequencies Q/2tt ~ 10 GHz [9], 
and also for moving vortex arrays driven by a large superimposed dc trans- 
port current (/ ^> Iq) [10]. At low frequencies (and 1 = 0) the linear ac 
response is strongly altered by pinning, so that a sample in the mixed state 
behaves as a true loss-free superconductor, with A" ~ [11] and a quasistatic 
penetration depth A' ~ 1-100 /im [12]. The way pinning and flux flow interfere 
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in the cross-over regime (Q/27T ~ 0.1-100 MHz in dirty conventional materi- 
als) is involved in the complex spectrum Z(Q), which consequently should be 
the main concern for the dynamics of pinning (see the theoretical review [13] 
and references therein). From a practical point of view it is equally impor- 
tant to know the width of the cross-over regime and the position of its center 
point, the depinning- frequency Q p , because they determine both the onset of 
losses at low frequencies and the recovery of the free flux-flow response at high 
frequencies. 

Our experiment is designed to study the skin-depth spectrum A ac (fi) in the 
cross-over regime. For this purpose, we measure the penetration of small rf 
magnetic fields into thick superconducting slabs in a large frequency range, 
n/2n ~ 2kHz-20MHz . 

The quasistatic regime is classically described in terms of an elastic inter- 
action between the vortex array and sample defects [14] (p. 370). Assuming 
the existence of a pinning force density, small vortex displacements u about 
their equilibrium position result from a Lorentz force J x w per unit volume, 
balanced by a linear restoring force F p = —a^u, where is the phenomeno- 
logical Labusch parameter. With increasing frequency a viscous term — tjvl 
(rj = (XfUJ 2 ) comes into play and damps vortex motion. Linearizing both the 
Josephson relation E = uj x v L (v L = ii) and the force-balance equation, one 
obtains a local constitutive equation in the form of the Ohm's law J± = a*E, 
where a* = o" f + ia L / (u 2 Q). This results in a one- mode electrodynamics with 
a complex penetration depth A ac given by 



1 1 2i 

— = -I/ZOO" = ~\2 72 



-ifi a*Q = - T2- (4) 



Here Ac = oj / (^oCXl) 1 ^ 2 is the Campbell penetration depth for small qua- 
sistatic vortex motions [12]. The spectrum (4) has been amended so as to 
include thermally-assisted vortex motion (flux creep) [15,16]. Thermal acti- 
vation entails an additional ^-dependence in «l (or Ac) at low frequencies, 
but it preserves the electrodynamic description in terms of a local effective 
conductivity a* [13]. 

According to a recent two-mode electrodynamics [17,18], which follows from 
linearized Eqs. (1) and (2), there is, apart from the ordinary flux-flow mode 
(FF-mode), a second mode associated with strain and corresponding to curie 
terms in Eq. (1): this vortex-strain mode (VS-mode) is evanescent and non 
dissipative. It is worth noting that in Eq. (1) (like for superfluids in [4]) shear 
effects are disregarded. We adopt a rather old line of reasoning [19], according 
to which the small shear rigidity of the vortex array may be ignored in most 
dynamic problems. It remains that the shear modulus is essential for ques- 
tions about the detailed equilibrium configuration of the vortex lattice; we 
return to this point in Sec. 6. On the other hand one cannot ignore the role of 
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the VS-mode: although it dies off over a short depth Ay > a from the surface, 
it can carry large non-dissipative currents, «J S = —curie (J s < e/\y ~ 10 7 
A/cm 2 ). On introducing the appropriate surface condition for small reversible 
vortex motion, the amplitude and the phase of the two modes, and ultimately 
the surface impedance of the sample, are determined. It turns out that for the 
ac response of a perfect sample the surface condition u\\ = applies, reducing 
the amplitude of the VS-mode to its (small) diamagnetic contribution [10]. 
However, as suggested recently [17,20], the VS-mode can be considerably en- 
hanced by surface roughness: surface irregularities on the scale of a allow for 
an hysteresis in the "contact angle" between the vortex array and the aver- 
aged smoothed sample surface (y\\ < 0.1). On substituting the vortex slippage 
condition vn = —u/l, for the ideal boundary condition, we find a number 
of important results. The hallmark of the two-mode electrodynamics is the 
original shape of the A ac (f2) spectrum (see Eq. (14) and Figs. 2 and 3 below). 
Furthermore, we predict unexpected size-effects which are discussed in Sec. 5 
and Fig. 6. 

In this paper we detail the experimental procedure (Sec. 2), show new results 
on Pbln alloys over an extended frequency range both for perpendicular and 
oblique fields, and discuss finite-size effects at some length. A preliminary 
investigation of the sheath-superconductivity regime (above Bq 2 )) included in 
Sec. 4, supports our interpretation of the mixed-state spectra in term of two 
modes. The paper ends with a discussion about pinning mechanisms. 



2 Experimental principle 

We have investigated a series of long superconducting slabs (length L in the 
x-direction, width w in the y-direction). The applied magnetic field B , in 
the 2/z-plane, makes an angle 9 with the ^-direction. We concentrate to the 
9 = and 9 = tt/4 orientations referred to below as normal and oblique 
field. A small ac ripple b e~ lUt along the x-axis is produced by a 10-100 turns 
excitation coil. Induced currents make the vortices oscillate in the xz planes. 
The ac flux <3>*e~* m through the cross-section w x t is detected by means of 
a 10-100 turns pick-up coil wound on the sample and measured with phase- 
sensitive lock-in analysers (EGG 5202 and 5302) in the broad Q/2ir =2kHz- 
20 MHz frequency range. In order to obtain a high resolution in the loss-angle 
ip = arctan(A"/A') (S(p < 0.3 degree below 1MHz) we used thin bronze wire 
to make pick-up coils transparent to the ac-field, and the rf chamber was 
enclosed in a high-conductivity copper can. At the low excitation levels of the 
linear regime (bo < 1 fiT in Fig. 1), the resolution was ultimately limited to 
5\ < 0.1 fim by electronic noise. Flux-flow conductivities and critical currents 
were measured beforehand from dc volt age- current characteristics (/ < 25 A). 
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Most measurements reported in this paper were carried out on Pblnl0.5%wt 
slabs 40 (L) x 6(w) x t mm 3 with t = 0.3, 1.3, 4 mm. Homogeneous lead-indium 
alloys are easily prepared from pure metals, and extensive measurements in 
the literature give reliable values of their B C2 , k and p n (10.3/Alcm) [21]. 
Slabs were spark-cut directly from ingots and optionally chemically etched. 
With k = 3.5 magnetization effects are negligible down to very low fields, 
and uj = Bq. On the other hand, Bc2 (0.29 T at 4.2 K) is low enough so that 
the field-phase diagram can be investigated up to the normal state with our 
0.8Tesla magnet (see Fig. 1). Measurements on very clean niobium (p n < 10 
nftem), first reported in Ref. [20], will be detailed elsewhere. We show in Fig. 3 
new data on Pbln5%u>t (p n = 6.5 pQcm) and standard vanadium (99.5% with 
p n = 0.6/ificm), which exhibit the same behaviour as Pblnl0.5%wt. 

For the sake of obtaining absolute values of A ac we need to calibrate the ac-field 
b accurately both in amplitude and phase. Several techniques can be used for 
this purpose: the phase can be obtained from that of the current in the rf 
excitation coil, or by measuring the magnetic flux either through an auxiliary 
pick-up coil, or through the gap between the pick-up coil and the sample 
in the Meissner state B = 0. The Meissner state is used throughout this 
work as the reference <3>* = state, neglecting the penetration over the small 
London depth Al ~ 0.1pm. The most accurate technique, however, makes 
use of the complex susceptibility of the normal state xl- Knowing p n , the 
flux distribution in an isotropic normal metal can be easily calculated thanks 
to the homogeneous boundary conditions achieved in this geometry. Note in 
this respect that the situation is far less tractable in resistivity experiments 
where the spatial distribution of the ac applied current is an intricate problem; 
moreover, the measured voltage strongly depends on the leads arrangement. 

For a precise calibration of the phase of the exciting field in the normal state, 
we need the exact two-dimensional solution for the field distribution in a 
rectangular cross section w x t, which reads : 

b 16 k 2 . miry . lixz 
— = } ^~To sm -sin 

u m,n odd ran I n 

2 m 2 iT 2 n 2 7r 2 

V k mn _,2 ^ p2 ' ' 



w 



Integration over the section yields the apparent susceptibility x n = ^n/^owt : 

64 k 2 

n m^ d mW k 2 mn -2i/5 2 - U 



The inset in the lower Fig. 1 shows x' n an d x'L as a function of 5 n /t oc l/Vil, 
calculated from Eq. (5) for w = 4t. As a cross check we always observed that 
the phase determined by other methods coincides with that obtained from the 



5 



normal-state response within 0.3 degrees below 1 MHz. The calibration from 
the normal-state response is, for a large part, unaffected by stray capacitances 
from which coil circuits usually suffer at high frequencies. 

Whereas a 2D calculation was necessary in the normal-state response, the 
situation is different in the mixed state because of the strong anisotropy of the 
critical currents and the flux flow resistivity. In particular, as can be checked 
experimentally by making 9 = ir/2, the flux penetration through the lateral 
faces parallel to the vortices is negligible below 0.9fi C 2- Flux penetration in 
normal field can thus be regarded as a one-dimensional problem with $* = 
2wbo\nc- Anyway, the question of the lateral penetration can be circumvented 
by working in oblique field where the 4 faces play symmetric roles so that 
$* = 2(w + t)b \ ac . In oblique field, bulk currents flow at an angle 6 = 7r/4 
with vortices, and the effective resistivity is p i (6) = p fJ _ cos 2 

The experimental procedure is the following: the flux $* = $(_B ) — $(0) is 
measured as a function of B at constant Q up to the normal state where the 
calibration is achieved. Then the data for each vortex state are collected as a 
function of to yield the complex spectrum. In Pblnl0.5%wt slabs the resolu- 
tion was not affected by the reproducibility of the vortex states. Otherwise we 
are careful to reproduce the same vortex state, either by using the same field 
history or by removing metastability by applying a large transient over-critical 
dc current. The frequency-range is limited to two decades for a given arrange- 
ment of exciting and pick-up coils; it is extended to four decades in Figs. 2 
and 3a by juxtaposing data taken from two distinct experimental setups with 
different windings, with several days lying between the two measurements. 



3 Skin effect in the mixed state 

Let us first recall the theoretical steps leading to the two-mode impedance 
spectrum (14), developed in Refs. [17,18,20]. 

Consider a superconducting half-space z < 0, immersed in a constant magnetic 
field B z, and submitted to a small transverse ripple boxe'^K In the absence 
of Hall effects (f3' = 0), vortex modes are linearly polarized, and the ac fields 
b x , E y , u x , u = v x /ik, . . . split into two modes, 




-int 



(6) 
(7) 




(9) 
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with respective wave numbers [18] 

l/ifci = A f = (1 + i)<J f /2, l/ifc 2 = A v = \y/(ioe/{u) + (ioe). (10) 



&i (resp. 62) and v\ (resp. v^) are the complex reduced amplitudes of the FF- 
mode (resp. VS-mode). Field continuity implies that 61 + 62 = 1- Additional 
relations between magnetic and vortex fields are obtained by linearizing the 
London equation (2) [18]: 

vi = h — ; u 2 = -b 2 — —. (11) 
-do Hq£ -do 

Note in Eqs. (11) the different behaviour of the vortices in the two modes: 
whereas in the FF-mode vortices follow the magnetic field (oui = B ui = 61), 
in the VS-mode vortices incurve in the opposite direction with a deflection 
drastically enhanced by the large factor u/fj, e ~ B /M >> 1, where M is the 
reversible magnetization of the superconductor. 

The apparent penetration depth of some combination of two modes is a 
weighted average of the individual penetration depths Af and Ay: 

A - = -^sr = ^iAf + My ^ 61 A f . (12) 

Since Ay % Al = A (a; = 0), the VS-mode can be regarded as superficial 
on the scale of the sample, so that its contribution to A ac in (12) can be 
neglected within our experimental resolution (~ 0.1 /im in A' or A"). The 
surface mode intervenes not only as a discontinuity in the field amplitude, but 
also, and perhaps less intuitively, as a phase-lag (arg&i) in the bulk mode. The 
latter effect turns out to be important in understanding the low-frequency ac 
penetration. It will be seen below that both amplitude and phase of b\ are 
controlled by the boundary conditions for vortices at the surface. 

Let us now introduce the boundary conditions for a rough surface and small 
vortex displacements (u < 1 A) [20]: 

^x(O) = v\ + v 2 = -wi/7. (13) 

It should be noted that Eq. (13) is independent of the frequency. This comes 
from the fact that any distortion of the vortex lattice matching the surface 
roughness in the small depth Ay can be regarded as quasistatic motion as 
far as <5f ^> Ay i.e. in the whole investigated range of frequencies. The surface 
roughness is characterized by the length I. A geometrical interpretation of I can 
again be borrowed from superfluids: an isolated vortex, pinned at the bottom 
of a spherical bump, when acted on in the bulk (e.g. by a second-sound wave), 
bends near the wall so as to keep on ending normal to the surface; whence a 
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linear relation v x = —u/l with I = R, the radius of curvature of the bump. 
Eq. (13) is nothing but a generalization of this picture to the collective motion 
of vortices interacting with their images in the disordered surface. 

By extrapolating Eq. (13) for large displacements u ~ a we obtain the maxi- 
mum standard deviation of vortices at the surface, v x ~ a /I. This value of v x 
also determines the surface critical-current density Kq = ev x ~ ea/l accord- 
ing to the Mathieu-Simon model of the critical state [1]. From experimental 
values of the critical current we infer that v x ~ 10 _2 -10 _1 , and I ~ 1-10 jum. 

Using condition (13) together with equations (8), (9) and (11), we determine 
b\ in (12), which gives the following expression for the effective penetration 
depth: 

1 11 — i . r r lu 
— = — + —=— , where L s = I (14) 

Aac -^S Of Hoe 

is a length directly related to the roughness length /, and therefore connected 
with the order of magnitude of critical currents by Kq ~ ua/ po-^s- 

The frequency dependence of A ac is governed by that of the flux-flow skin 
depth 5{ oc 1/y/U. The right-hand side of Eq. (14) is a straight line in the 
complex plane, which, by inversion, gives a quarter-circle for the spectrum 
^ac(^) (Fig- 3). Through the action of the boundary conditions, it turns out 
that, in the low- frequency limit, A ac ~ 6iAf ~ L s is real and independent of 
fl The length L s , as A'(0), plays the role of the Campbell length. We stress 
however that, in contrast to Ac, Ls is not an actual penetration depth, but 
rather some fraction of a much longer depth 5{. 

Fig. 2 shows a typical complex spectrum of A ac in a log-frequency scale. Com- 
pared to our previous measurements [20], the number of frequency points and 
the frequency span have been doubled so as to cover most of the depinning 
transition. We confirm our observation that the cross-over is much broader 
than predicted by the Campbell bulk-pinning expression (4) (dotted lines in 
the figure) and extends over 6 decades of frequency (10 3 -10 9 Hz in Fig. 2). This 
observation is in apparent contradiction with that of Gittleman and Rosen- 
blum [11] who reported a depinning transition in the ac resistivity of similar 
materials, over 2-decades of frequency. In fact, there is no inconsistency if 
one takes into account the finite thickness of their samples (see Sec. 5 be- 
low). The quality of the fit to Eq. (14) (solid line) in Fig. 2 is, in itself, a 
strong indication for the two-mode electrodynamics. Apart from theoretical 
considerations, formula (14) proves powerful in extracting reliable values of the 
free-flow resistivity pf within a few percent. This method compares favorably 
with microwave techniques. 

The universal shape of the depinning spectra is well illustrated in the Argand 
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diagram of Fig. 3, where data from two Pbln alloys (p n = 6.5 and 10.3 pficm 
respectively)) and a vanadium sample (p n = 0.6 /jficm) have been collected. 
The characteristic quarter-circle is observed in a variety of experimental con- 
ditions and samples, including very pure niobium (p n < lOnficm) in Ref. [20]. 
Two distinctive features, well observed in Fig. 3 a, are worth mentioning : i) 
the linear behavior A' + A" = A'(0) at low frequency, and ii) the maximum of 
A", A" = |tan(7r/8)A'(0) = 0.207A'(0), which is 30% less than the generally 
accepted classical value 0.272A'(0). 

The field dependence of both L s and critical current densities K c are shown 
in Fig. 4 (upper). The order of magnitude of the roughness length / can be 
obtained by using the Abrikosov expression fi^e = (Bc2~ uj) /2.32k 2 in l/L$ = 
Hqe/Iuj and u = B . With n = 3.5 we find reasonable values 0.4 < I < 1 jum 
(dashed lines in Fig. 4). Moreover, K c = ic/2u> values, deduced from dc 
measurements, are also found to agree with the estimate K c ua/ fi L s . Fig. 4 
(lower) illustrates the resolution in pf(u>) achieved with rf surface impedance. 



4 Skin effect through the superconducting sheath 

We have investigated the field range B /Bq2 = 0.9-1.1 to explore the be- 
haviour of the ac response in the transition from the bulk mixed state to the 
surface superconductivity. As surface superconductivity is not destroyed in 
oblique field, we took advantage of this geometry to symmetrize the role of 
the four faces of the sample (see Sec. 2). We recall that the superconducting 
sheath (depth ~ £) is populated by short vortices (flux spots) with a density 
u = B j_ (see Ref. [22] and references therein). 

The most striking fact in Fig. 5 is that, at the transition, there is no change in 
the shape of spectra, which, both above and below Bc2, obey the theoretical 
expression (14). Above Bq 2 the bulk is normal, and the existence of screening 
surface currents is obvious; therefore, the two-mode electrodynamics and the 
free bulk penetration of the wave, which underlie Eq. (14), are not surprising. 

Concerning the field-dependent parameters involved in Eq. (14), Ls and pf, we 
observe no discontinuity of Ls(B ), whereas a jump is observed in pf as shown 
in Fig. 5b. The continuity of Ls(B ) is consistent with the observation that 
critical currents decrease monotonically through the transition together with 
the strength of the surface vortex state. The vanishing of the vortex structure 
in the bulk at Bq 2 results in an abrupt jump from p n cos 2 6* = p n /2 to p n as 
expected from the anisotropy of resistivity in the mixed state [22] . 
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5 Size effects in thin slabs 



Equation (14), and results in the above Section, concerned the ac response 
of two independent half spaces. In thin slabs and/or at low frequencies, fields 
penetrating from opposite faces interfere in the bulk whenever t < 2Sf . Assum- 
ing symmetric surface conditions on both faces z = ±i/2, the corresponding 
one-dimensional solution takes the even form : 



b ~ cosh(z/Af) 
k) ~ 1 cosh(t/2A f ) 



+ ~b 2 cosh ^/ Av ) 



E 



-iQb 
bo 

60 



- iAf sinh(,/A f ) ^ 



cosh(t/2A v ) 
sinh(z/Ay) 



cosh(t/2A f ) 

~ cosh(2;/Af 
= b\ 



cosh(t/2A v ) 

uj ~ cosh^/Ay) 
62- 



Mi 



cosh(t/2Af) jj, e cosh(t/2Av) 
sinh(^/Af) uj ~ sinh(z/Av) 
cosh(t/2Af) jj e 2 V cosh(t/2Av) 



(15) 
(16) 
(17) 
(18) 



The effective penetration depth reads 

A ac = E ^ Z^y^^ = Mf tanh(t/2A f ) + 6 2 A V ^ M f tanh(t/2A f ). (19) 
— 2iilb 

Thin films (t < 2Ay) are excluded here by assuming tanh(t/2Ay) = 1. Inserting 
boundary conditions v{±t/2) = ^fu 1 (±t/2)/l, Eq. (14) becomes 



1 


1 




1 
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r**"^ (20) 

(t/5t -> 0). (21) 



Comparing Eq. (20) with Eq. (14), it is seen that size effects will start when t 
is decreased below about 35 f , where the coth factor deviates significantly from 
unity. As a result losses are significantly reduced, the loss angle (p falling off 
below the half-space level given by siny?= |A ac |/5f- Data (</?, |A ac |) obtained at 
100 kHz in Pbln slabs, where Pf(B ) (then 5{) had been separately measured, 
clearly demonstrate these size effects (see Fig. 6). The fact that for decreasing 
sample thicknesses size effects are encountered when t ~ 5f, i.e. well before 
t ~ L s , is not consistent with any one- mode theory (see Sec. 1). 

According to Eq. (21), a second size effect occurs in the thin limit when 
t ~ 2L S . This second "dimensional" cross-over, which affects now the real 
part A', could wrongly suggest that L$ is a true length in accordance with 
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one-mode approaches. This shows that one must be careful in interpreting 
susceptibility measurements. 



Let us conclude by briefly examining resistivity measurements in foils within 
the thin limit. We refer in this respect to the pioneering work of Gittleman 
and Rosenblum [11], as also to recent measurements in YBaCuO platelets [23]. 
The analysis of such resistivity spectra is much more intricate. Many authors, 
using this method assume that the current density does not vary along the 
width w of the foil; this requires the absence of two-dimensional skin-effects, 
S{ 3> w > t, a condition that is hardly fulfilled over several frequency decades. 

The one-dimensional treatment of the ac resistivity relies on the odd solution 
of Eqs. (6)-(9). The effective resistivity of the foil is defined as 

£(V 2 ) t . ^ t 

Pac = Po — 7 o = -w2A ac -, (22) 

o 2 2 



where 



— = — + -— tanh^— -L. (23) 
A ac L s d { 2d { 



In the limit t < 25 u Eqs. (22), (23) reduce to : 

-ifiott 2 2 2% 

Pac £A ac tLs 8 { 



2 ' 



(24) 



Surprisingly enough, we find the same resistivity spectrum as that given 
by Eq. (4) which follows from the bulk-pinning theory. Both Eqs. (4) and 

(24) equally fit the Gittleman-Rosenblum data (by setting Ac = ^jL s t/2 in 
Eq. (24)). We are led to the conclusion that resistivity experiments alone can- 
not decide between very contrasted models of pinning. 



6 Conclusions 



We have reported the first comprehensive investigation of dynamic depinning 
transition in the surface impedance geometry. It relies on the description of the 
depinning transition through Eq. (14) and requires the use of thick samples 
(t > 35f). The agreement is so good that surface-impedance measurements in 
the rf-domain provide a precise and accurate way of measuring both parame- 
ters Ls and pf. As shown in Fig. 4, dc data on the flux-flow resistivity, when 
available, confirm the accuracy of the method; in other cases where low-level 
ac techniques are necessary (large critical currents, ultra-low temperatures, 
. . . ), the surface impedance becomes the best suited non-invasive method for 
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measuring pf . Apart from this metrological aspect, it is clear that, as concerns 
the pinning length L$ = A'(0), our interpretation is completely at variance 
with classical views on pinning ([14] p. 345). In this respect we wish to add 
some comments. 

The remarkable agreement with the two-mode spectrum (14) as well as the 
observation of predicted size effects for t ~ 5{ at low frequencies bring out 
two important facts: i) The effective length, as deduced from the surface in- 
ductance X/VL = /i A'(0), is not an actual penetration depth; ii) free flow is 
observed in the bulk (over thicknesses t > 1 mm), in spite of a large density 
of volume defects in a variety of polycrystalline samples from alloys to pure 
metals. 

In the literature, A'(0), usually referred to as the Campbell penetration depth, 
is directly expressed in terms of the bulk-pinning Labusch parameter a L (see 
Sec. 1). In soft samples, such as used in this work, bulk pinning is commonly 
analysed in the frame of the collective pinning theory of Larkin-Ovchinnikov 
[14], where is in turn related to Larkin transverse and longitudinal cor- 
relation lengths i?c and L c . These correlations lengths are intricate expres- 
sions of thermodynamic and pinning parameters (elastic constants c 44 and 
C66, core radius £, pinning-centers concentration n and the individual pin- 
ning force /). For thick samples (3D collective pinning) one expects that 
A'(0) ~ Lc ~ C44C66C 2 /nf 2 [13]. Thus many authors consider the lengths 
A'(0) deduced from the low-frequency susceptibility as a routine measurement 
of L c and merely investigate the field and temperature dependence of L c . 
This studies must be taken with due caution; our results show that a spectral 
analysis should precede any interpretation. 

The role of point defects and elastic constants is certainly important in ana- 
lyzing the short and long-range distortions of the vortex lattice, as observed 
in a number of techniques (Bitter decoration, neutron scattering, STM, ... ). 
However, our experiments seriously question the relevance of this mechanism 
of weak pinning in transport properties of soft samples. Here, we indeed ex- 
clude the case of hard materials where additional interfaces in the bulk may, 
to some extent, act like the surfaces of soft samples (twin planes in YBCO 
[24] , precipitates in industrial wires, ... ) . 

The existence of a free vortex flow in the bulk, and the exclusive role of the 
surface in the critical-current properties of soft superconductors, have been 
confirmed in previous experiments under over-critical conditions: for instance 
the voltage and field noise in the flux- flow regime [25], the high-level step- 
response of a slab in parallel field [26], or surface localization of sub-critical 
currents [27], . . . We emphasize that these experiments have been specially 
designed to discriminate surface and bulk currents. It is clear that the existence 
of bulk sub-critical currents would imply bulk pinning, but there is up to now 
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no experimental evidence for such bulk currents in soft samples. In particular 
one cannot, in this respect, rely on the recent Hall probe technique, which 
determines the radial distribution of the current in platelets. The observed 
pile-of-sand profile of the magnetic field confirms the critical state model [28], 
but critical state model should not be confused with bulk-pinning model. The 
critical-current densities in this geometry may be superficial and the technique 
used cannot resolve the axial distribution. 
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Fig. 1. Transverse susceptibility of a 6.1(u>) x 1.22(t) mm 2 Pblnl0.5%wt. slab 
(p n = 10.3 ^$7cm) in a broad range of normal fields including the mixed state 
(B < 0.29 T) and the surface sheath regime (0.29 T< B < 0.49 T). The upper 
panel shows the quasistatic response at low frequency with vanishing losses (%") 
and a finite in-phase apparent penetration depth A' = x't/2. In the lower panel, Q 
is close to the dep inning threshold Q p , and \" ls clearly visible. For each frequency, 
the normal state response (plateau for Bq > 0.49 T) is used to calibrate the phase 
and the amplitude of the excitation by comparison with a two-dimensional calcula- 
tion (inset of the lower panel and Eqs. (5)). Non-linear effects are displayed in the 
inset of the upper figure: strictly the linear response requires extrapolating data to 
bo = for each field and frequency; in practice we work with a constant but small 
level bo ~ 1 /uT (or induced currents K = bo/ Ho ~ lOmA/cm), which translates into 
minute vortex displacements u ~ X'bo/Bo ~ 1 A. 
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Fig. 2. Broad-band spectrum of the apparent penetration depth A ac in a 
6(w) x 4(£) mm 2 Pblnl0.5%u> slab in oblique field. It shows the cross-over be- 
tween the quasistatic regime (Q/2ir < 1 kHz) and the free flux-flow regime 
(Q/2ir > 1 GHz). The 2 kHz-20 MHz experimental range has been achieved by using 
two different set-ups, each covering two frequency decades. Solid lines are theoreti- 
cal fits to Eq. (14) with Ls = 42 ± 0.5 fim and p{ = 4.0 ±0.1 fiQcm. For comparison 
the Campbell spectrum (4) (dashed- line), with the same set of parameters, yields a 
much narrower transition with a maximum A" /Ls = 0.272 larger than the maximum 
\"/L s = 0.207 of the solid line. 
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Fig. 3. The apparent penetration depth A ac (fi) displayed in the complex plane. Data 
have been measured in the mixed state of conventional isotropic superconductors: 
Arc of circles (solid lines) are fitted to the data according to Eq. (14) by adjusting 
the skin length L$ (real-axis intercept indicated by an arrow in (a)). A typical bulk 
pinning spectrum, given by the Campbell expression (4) with Ac = 63 //m, has been 
plotted (dotted line) for comparison to the 0.24T data, (a) a 6(w) x 4(t)mm 2 slab 
of Pbln(10.5%iui) (p n = 10.3/iO.cm); (b) a 5.8{w) x 2(t)mm 2 slab of Pbln(5%wt) 
(Pn = 6.5/ifi.cm); (c) a 6.7(w) x l(t)mm 2 vanadium slab (p n = 0.6/iJ7.cm), 
spark-cut from a cold-rolled commercial plate. 
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Fig. 4. Upper panel : Typical field-dependences for the skin-length Ls (bullets) 
and the surface critical current density Kq (squares) in Pblnxo^u, (-B C 2 = 0.29 T, 
k ~ 3.5) at T = 4.2 K. Dotted lines are theoretical curves Ls = IBq/hqe us- 
ing the Abrikosov expression e ~ (B C 2 — Bq) /2.32/Uok 2 , and a pinning length 
I = 1.1 /um (resp. / = 0.4 /xm) for intermediate fields (resp. high fields). The prod- 
uct KqLs ~ 10 mA agrees well with the theoretical estimate ^/Bipo/ /xo (= 14 mA 
at So = 0.15 T) for surface pinning. Lower panel : the field dependence of the 
flux flow resistivity in Pblnl0.5%u> at T = 4.2 K deduced from the slope of dc 
voltage-current characteristics (squares) and from rf surface impedance (bullets), rf 
data are deduced from the frequency dependence of A ac (f2) in Fig. 3a. 
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Fig. 5. Upper panel: Continuity of the spectra A ac (S7) across the 
bulk-superconductivity transition in oblique field. Data are taken from a 
6.1(u>) x 1.22(£) mm 2 Pblnl0.5%u>t at T = 4.2 K (,602 = 0.29 T). In this geometry, 
currents flow at 45 degrees to the magnetic field, and the penetration of the ac-field 
is symmetric over the four faces of the slab throughout the investigated range 
0.9 < Bq/Bq2 < 11- Agreement with the theoretical spectrum (14) is observed 
on both sides of the transition which is a strong indication that the screening 
mechanism is the same in both cases. Lower panel : The skin-length L${Bq) (•) 
and the flux-flow resistivity Pf(Bo) (o), deduced from the fits of A ac . Ls(Bo) is itself 
continuous through the transition, as is the dc critical current. The only accident 
observed at Bq2 is a drop by a factor 2 in p^ir/A) which signals the vanishing of 
the mixed-state anisotropy Pf(6) = p{ (0) cos 2 0. 
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Fig. 6. The sine of the loss- angle sin 99 = A"/|A ac | as a function of |A ac |/(5f in two 
slabs of the same Pblnl0.5%ti;i material with different thickness t = 0.3 mm (•) and 
t = 1.3mm (o). The skin-depth S{, or the flux-flow resistivity pf, is deduced from 
separate dc voltage-current characteristics. Data points (|A ac |,(/j) have been taken 
at a constant low frequency as a function of magnetic field : Bq = 0.04-0.18 T and 
S t = 0.13-0.32 mm in the t = 1.3mm slab; B = 0.06-0.27T and S t = 0.15-0.45 mm 
in the t = 0.3 mm slab. The theoretical law sin^ = |A ac |/5f for the half space is well 
obeyed in the thick slab with Sf < t/A, while a significant deficit in the losses of thin 
slab is observed for 5{ ~ 0.66-2t. Note that, |A ac | ~ 5-30 /im and t = 10-60|A ac | in 
the same field range, so that no size effect would be expected if A ac was the actual 
penetration depth. Quantitative comparison with theory is difficult in this thin 
limit as losses become drastically sensitive to the symmetry of surface conditions 
(compare Eq. (19) with Eq. (23)). 
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